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Abstract 

A  mathematical  model  to  calculate  the  volume  specific  three-phase  boundary  length  in  the  porous  composite  electrodes  of  solid-oxide  fuel  cell 
is  presented.  The  model  is  exclusively  based  on  geometrical  considerations  accounting  for  porosity,  particle  diameter,  particle  size  distribution, 
and  solids  phase  distribution.  Results  are  presented  for  uniform  particle  size  distribution  as  well  as  for  non-uniform  particle  size  distribution. 

©  2007  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Three-phase  boundary  (TPB)  is  a  geometrical  parameter  that 
is  of  crucial  importance  for  the  performance  of  solid-oxide  fuel 
cells  (SOFC)  [1].  In  general  a  high  TPB  is  required  for  high 
electrochemical  performance.  Even  though  there  are  number  of 
publications  concerning  the  modeling  of  SOFCs,  most  of  these 
models  lump  the  TPB  length  in  the  exchange  current  density 
parameter,  which  makes  the  modeling  task  much  easier  [2-5]. 
Nevertheless,  it  is  important  to  understand  the  influence  of  elec¬ 
trode  micro- structure  on  the  TPB  length.  Deng  et  al.  [6]  reported 
a  mathematical  model  to  calculate  the  TPB  length  for  mono¬ 
sized  particles.  Their  expression  leads  to  a  TPB  length  with  units 
of  m/m4.  It  is  quite  difficult  to  conceptualize  a  physical  quantity 
with  the  above  mentioned  units.  Furthermore,  their  model  do 
not  account  for  the  coordination  number  between  the  ionic  and 
electronic  conductors.  In  most  cases  of  a  composite  electrode, 
one  is  interested  in  the  volume- specific  TPB  length. 

In  an  interesting  article  Wilson  et  al.  [7]  reported  the  three- 
dimensional  reconstruction  of  a  SOFC  anode  using  ion-beam 
scanning  electron  microscopy  to  determine  the  micro- structural 
properties.  They  estimated  the  volume-specific  TPB  length  for 
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their  sample  to  be  4.28  x  1012  m/m3.  Brown  et  al.  [8]  reported 
that  Ni  constituted  the  fraction  of  larger  particles  with  the  par¬ 
ticle  sizes  raging  0.5-3.0  p,m,  while  YSZ  phase  particle  size 
distribution  were  0. 5-1.0  pan.  So  it  is  also  important  to  consider 
the  volume  specific  TPB  length  for  cases  with  two  different 
mono-sized  particle  distributions.  Schnieder  et  al.  [9]  reported 
the  discrete  modeling  of  composite  electrode  and  developed  an 
analytical  model  for  calculating  the  TPB  length.  They  reported 
that  the  maximum  TPB  length  is  obtained  when  the  volume 
fraction  of  ionic  conducting  particles  is  50%.  They  also  devel¬ 
oped  an  analytical  model  for  calculating  the  TPB  length,  which, 
however,  predicts  non-zero  coordination  number,  even  when 
the  ionic  phase  volume  fraction  is  zero.  Similar  to  the  anode, 
the  three-phase  boundary  length  in  cathode  is  important  for  the 
electrochemical  reduction  of  oxygen.  Chan  et  al.  [10]  reported 
an  analytical  expression  for  the  calculation  of  TPB  length.  In 
this  work  we  present  the  development  of  a  mathematical  model 
to  describe  the  volume-specific  TPB  length  per  unit  electrode 
volume  from  a  geometrical  perspective. 

2.  Model  development 

The  model  presented  here  is  developed  by  considering  a  geo¬ 
metric  volume  of  a  composite  electrode  characterized  by  its 
porosity  0  and  particle  radii  r\  and  r2.  It  is  assumed  that  the  ionic 
and  electronic  solid  phases  are  made  up  of  spherical  particles. 
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Nomenclature 

a 

radius  of  lens  (m) 

d 

distance  between  the  sphere  centers  (m) 

^tpb 

average  TPB  length  (m) 

M 

ratio  of  the  number  of  ionic  and  electronic  parti¬ 
cles 

Nv 

total  number  of  particles 

% 

number  of  particles  with  radii  n 

Np2 

number  of  particles  with  radii  r2 

r 

radius  (m) 

F^tpb 

volume  specific  TPB  length  (m/m3) 

Vl 

volume  loss  (m3) 

vt 

total  volume  (m3) 

Vti 

total  volume  loss  (m3) 

z 

average  coordination  number 

Ze 

coordination  number  of  electronic  conductor 

Z; 

coordination  number  of  ionic  conductor 

Zi— e 

coordination  number  between  ionic  and  elec¬ 
tronic  conductor 

Greek  letters 

a 

mean  radius  ratio 

cp 

gas-phase  porosity 

0e 

volume  fraction  of  electronic  phase 

0i 

volume  fraction  of  ionic  phase 

xj/ 

fractional  overlap 

Subscripts 

e 

electronic 

i 

ionic 

Here  r\  is  the  radius  of  the  larger  sphere,  r 2  that  of  the  smaller 
one,  and  d  is  the  distance  between  the  centers  of  two  spheres.  (In 
the  analysis  presented  here  we  assume  the  electronic  conductors 
are  larger  than  the  ionic  conductors  [8] .)  Further,  the  intersection 
of  the  spheres  is  a  curve  lying  in  the  plane  parallel  to  the  x-y 
plane,  whose  radius  a  is  given  by 


The  intersection  of  two  spherical  particles  is  only  possible  under 
the  following  constraint  which  is  derived  from  Eq.  (2) 


d  <  r\  +  r2.  (3) 

In  a  porous  electrode  several  combinations  of  the  spherical 
particles  are  possible.  For  the  sake  of  clarity  few  examples  are 
considered  here.  In  Fig.  2(a)  there  are  10  spherical  particles  and 
1 1  intersection  volumes,  and  each  intersection  leads  to  different 
volume  losses  depending  on  d.  If  the  value  of  d  entering  in  Eq. 
(1)  is  the  average  value  for  the  porous  media,  then  for  the  case 
shown  in  Fig.  2(a),  the  total  volume  loss  is  11  x  V\\  where  Vi 
is  the  average  volume  loss.  For  the  case  in  Fig.  2(b),  there  are 
eight  particles  and  eight  intersection  volumes.  However,  some 
intersection  volumes  are  overlapped  by  the  adjacent  intersection, 
in  this  case  the  total  volume  loss  Vti  is  expresses  as 

Vti  =  Np(l  -  f)V 1,  (4) 

where  x//  is  the  fractional  overlap  and  Np  is  the  number  of  par¬ 
ticles.  When  x/r  =  0,  there  is  no  overlap,  and  the  case  becomes 
similar  to  that  discussed  in  Fig.  2(a).  There  may  also  be  cases 
where 

Vti  =  (Np  -  1)(1  -  ir)V 1,  (5) 


We  start  the  model  development  by  considering  the  intersection 
of  two  spherical  particles  as  shown  in  Fig.  1.  The  volume  of  the 
three-dimensional  lens  common  to  both  the  spheres  as  a  result 
of  intersection  is  given  by 


7x{r\  +  r2  —  d)2(d2  +  2  dr2  —  3  r2  +  2dr\  +  6r2ri  —  3  r2) 
12 d 


(1) 


Fig.  1.  Intersection  of  two  spherical  particles. 


Fig.  2.  Two  possible  combinations  of  spherical  particles  in  a  composite  elec¬ 
trode. 
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as  in  the  case  of  a  three  sphere  intersection  with  two  intersection 
volumes.  However,  for  a  homogenized  porous  media  we  con¬ 
sider  the  case  given  by  Eq.  (4).  Let  Vt  be  the  total  volume  of  the 
electrode  under  consideration,  then 


particles  surrounding  an  electronic  particles  and  vice  versa  can 
be  written  as 

Zi_ e  =  Ze_i  =  0i0eZiZe.  (14) 


(1  -  m 't  =  Npl  i rcrl  +  Npf-7tr\  -  (Npl  +  Np2)(l  ~  f)Vx , 


(6) 


where  0  is  the  gas-phase  porosity,  Api  the  number  of  particles 
with  radius  r\ ,  and  Ap2  is  the  number  of  particles  with  radius  r2. 
If  M  is  defined  as  the  ratio  between  the  number  of  particles 


M=Uv2  =  <h± 

NPi  4>i  if 


(7) 


where  0i  and  02  are  the  volume  fraction  of  particles  with  radii 
r\  and  7*2,  respectively.  Then  Api  can  be  calculated  from  Eq.  (6) 
as 


Nvi  = 


_ (1  ~  0)Vi _ 

(4/3)7r(r3  +  Mr\)  -  (1  -  ^)(1  +  M)V i ' 


The  total  number  of  particles  Nv  is  then  given  by 


(8) 


The  total  TPB  length  per  unit  volume  or  the  volume- specific 
TPB  r’Ltpb  length  is  then  given  by 


F^tpb  — 


0(1  -  0)(1  +  M)Zi_e27ra 

(4/3):r  (r3  +  Mr3)  -  (1  -  ifl  +  M)V \ 


(15) 


In  the  case  of  mono-sized  particles  with  equal  volume  distri¬ 
bution  of  ionic  and  electronic  phases  the  total  number  of  particles 
is  given  by 


Vp  —  iVpi  +  Np2  — 


(i  -  <m 

(4/3)tt r3  -  (1  -  ir)V\ 


(16) 


The  volume-specihc  TPB  length  (vLtp b)  is  then  given  by 


_  <P(  1  -  0)Z;_e2 na 

tpb  “  (4/3):rr3  -  (1  -  f)Vx 


(17) 


3.  Results  and  discussion 


Np  —  Npi  -T  Np2  — 


_ (l-0)Vt(l  +  M) _ 

(4/3)tt (r3  +  Mr\)  -  (1  -  ^)(1  +  M)V\ 

(9) 


The  length  of  boundary  due  to  the  intersection  of  two  spheres  is 
given  by  the  circumference  of  a  circle.  This  becomes  a  “three- 
phase”  boundary  length  only  if  the  intersecting  particles  are  of 
different  phases  and  the  intersection  is  associated  with  a  pore 
space.  Therefore,  the  average  TPB  length  of  a  composite  elec¬ 
trode  is  give  by 


^tpb  —  A^pZi_e027Tf2, 


(10) 


where  Zi_e  is  the  co-ordination  number  between  the  ionic  and 
electronic  conductors.  That  is  the  number  of  ionic  particles  sur¬ 
rounding  a  given  electronic  particle.  According  to  Bouvard  and 
Lange  [1 1]  the  coordination  number  of  electronic  conductor  Ze 
in  a  mixture  of  ionic  and  electronic  conductors  is  given  by 


Ze  =  3  + 


Z  —  3 

0e  +  0i«2  ’ 


(ID 


and  that  of  the  ionic  conductor  is  given  by 


Z;  =  3  + 


(Z  -  3)a2 
<Pe  +  0i«2  ’ 


(12) 


Here  Z  is  the  average  coordination  number,  0i  the  volume  frac¬ 
tion  of  the  ionic  particles,  0e  the  volume  fraction  of  the  electronic 
particles,  and  a  is  the  mean  radius  ratio  of  ionic  to  electronic 
conducting  particles,  i.e. 


a  — 


n 


(13) 


From  simple  statistics  the  average  coordination  number 
between  these  two  kinds  of  particle;  i.e.  the  number  of  ionic 


3.1.  Uniform  particle  size  distribution 

The  following  results  are  presented  for  mono  sized  particles 
for  both  ionic  as  well  as  for  electronic  phases.  Quite  obviously 
Eq.  (17)  predicts  maximum  volume  specific  TPB  length  (uLtpb) 
at  50%  porosity  while  other  parameters  are  fixed.  Fig.  3  displays 
the  influence  of  grain  size  on  r>Ltpb  as  a  function  of  gas-phase 
porosity  0  for  equal  volume  distribution  of  ionic  and  electronic 
particles.  As  predicted  by  Eq.  (17),  maximum  vLipb  is  observed 
at  50%  porosity  for  all  the  grain  sizes  considered  and  the  r>Ltp b 
increases  with  decreasing  grain  size.  For  the  results  presented  in 
Fig.  3,  the  average  distance  between  the  particles  is  assumed  to 
be  90%  of  the  particle  radius  and  a  coordination  number  (Zi_e) 
of  7.2  is  used.  From  the  order  of  magnitude  the  results  are  in 
good  agreement  with  the  experimental  evaluation  of  vL tpb  by 
Wilson  et  al.  [7]. 


Fig.  3.  Volumetric  TPB  length  as  a  function  of  gas-phase  porosity  for  various 
grain  sizes  for  mono-sized  particles  with  equal  volume  distribution  (cf.  Eq.  (17)). 
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Fig.  4  displays  the  pLtp b  for  various  porosities  as  a  function 
of  grain  size.  For  these  calculations  the  distance  between  the 
particles  is  kept  at  90%  of  the  particle  radius.  Like  in  the  previ¬ 
ous  case  maximum  pLtp b  is  for  smaller  grains  and  increase  with 
increasing  porosity.  However,  it  should  be  noticed  that,  increas¬ 
ing  the  porosity  above  50%  will  bring  down  the  r»Ltp b,  as  it  is 
obvious  from  Eq.  (17). 

The  influence  of  the  average  distance  between  the  particles 
on  vLtp b  is  shown  in  Fig.  5.  For  the  calculations,  the  grain  size 
is  kept  constant  at  2.5  pm.  As  the  distance  between  the  particle 
increases  the  r>Ltp b  decreases.  Quite  obviously  when  the  distance 
between  the  particle  increases  the  radius  of  the  plane  formed  by 
the  intersection  of  the  two  spheres  decreases,  which  leads  to 
shorter  length  for  the  curve  and  hence  shorter  r>Ltp t>.  As  evident 
from  Eq.  (17).  the  r>Ltp b  increases  linearly  with  coordination 
number,  which  is  shown  in  Fig.  6.  Fig.  7  displays  the  influence 
of  over  fraction  if/  on  pLtp t>.  As  expected  the  pLtpb  is  high  for 
low  over  fractions. 

3.2.  Non-uniform  particle  size  distribution 

Fig.  8  depicts  the  vL tpb  against  the  volume  fraction  of  the 
ionic  phase  (<p\)  for  various  particle  diameters  for  the  ionic 


Fig.  6.  Volumetric  TPB  length  as  a  function  of  coordination  number  for  a  grain 
size  of  2.5  |xm. 


phase.  The  particle  size  of  the  electronic  phase  is  kept  con¬ 
stant  at  3.0  p,m,  and  an  average  coordination  number  (Z)  of  6 
is  assumed.  Generally,  the  pLtp b  increases  with  decreasing  par¬ 
ticle  diameter.  However,  the  porosity  for  the  maximum  pLtp b 


Fig.  5.  Volumetric  TPB  length  as  a  function  of  the  distance  between  the  particle  Fig.  8.  Volumetric  TPB  length  as  a  function  of  the  volume  fraction  of  ionic 
centers  for  various  porosities.  conductor  for  various  grain  sizes  of  the  ionic  conductor. 
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Fig.  9.  Total  solid  phase  volume  as  a  function  of  the  volume  fraction  of  ionic 
conductor  for  various  grain  sizes  of  the  ionic  conductor. 


Fig.  10.  TPB  length  as  a  function  of  the  volume  fraction  of  ionic  conductor  for 
various  grain  sizes  of  the  ionic  conductor. 

is  different  for  different  volume  fractions.  It  may  be  interesting 
to  note  that  the  vLtp b  is  higher  for  higher  grain  size  at  lower 
volumetric  composition  of  ionic  phase.  This  is  mainly  because 
at  lower  ion  conductor  volume  fractions,  the  total  volume  of  the 
solid  phase  is  lower  with  higher  grain  size,  while  the  three  phase 
area  is  nearly  the  same.  The  total  solid  phase  volume  of  higher 
grain  size  is  lower  because  of  the  higher  volume  loss  (Vi).  The 
total  solid  phase  volume  as  a  function  of  volumetric  composi¬ 
tion  is  displayed  in  Fig.  9,  and  the  three-phase  area  is  shown  in 
Fig.  10. 

Fig.  1 1  displays  the  vLtp b  for  two  different  cases,  a  case  where 
the  particle  size  of  the  ionic  phase  is  larger  than  that  of  the  elec¬ 
tronic  phase,  and  a  case  where  the  particle  size  of  the  electronic 
phase  is  larger  than  that  of  the  ionic  phase.  When  the  electronic 
phase  particle  diameter  is  smaller  compared  to  the  ionic  phase, 
the  maximum  in  vLipb  occurs  at  lower  volume  fraction  of  the 
ionic  phase.  And,  when  the  electronic  phase  particle  diameter 


length. 

is  larger  compared  to  the  ionic  phase,  the  maximum  in  vLtp b 
occurs  at  higher  volume  fraction  of  the  ionic  phase. 

4.  Conclusions 

We  have  developed  a  mathematical  model  for  volume- 
specific  TPB  length  vLipb  based  simple  geometrical  concepts. 
The  model  takes  into  account  of  the  micro- structural  properties 
porosity,  particle  diameters,  and  volume  fractions  of  ionic  and 
electronic  phase.  The  geometrical  model  gives  result  which  are 
in  reasonable  agreement  with  experimental  evaluation  of  vLtp b 
[7].  The  model  presented  can  be  used  for  the  calculation  of  vol¬ 
ume  specific  TPB  length  for  uniform  particle  size  distribution 
as  well  as  for  non-uniform  particle  size  distribution. 
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